function [] = hyper(hint, nnit, npro)
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here

      if(nnit>0) 
            %     CALCULATIONS involving NITROGEN HYPERFINE
            if(nnit==1) 
                % Splitting by one nitrogen
                hint(1) = 0.333;
                hint(2) = 0.333;
                hint(3) = 0.333;
            elseif(nnit==2)
                % Splitting by two nitrogens
                hint(1) = 0.111;
                hint(2) = 0.222;
                hint(3) = 0.333;
                hint(4) = 0.222;
                hint(5) = 0.111;
            elseif(nnit==3)
                % Splitting by three nitrogens
                hint(1) = 0.037;
                hint(2) = 0.111;
                hint(3) = 0.222;
                hint(4) = 0.259;
                hint(5) = 0.222;
                hint(6) = 0.111;
                hint(7) =0.037;
            elseif(nnit==4)
                % Splitting by four nitrogens
                hint(1) = 0.012 ;                                      
                hint(2) = 0.049  ;                                     
                hint(3) = 0.123;
                hint(4) = 0.198;
                hint(5) = 0.235;
                hint(6) = 0.198;
                hint(7) = 0.123;
                hint(8) = 0.049;
                hint(9) = 0.012;
            end
      else
         %   couplings to protons or other spin-1/2
         
         %  splitting by 1 proton
        if(npro==1)
            hint(1) = 0.5;
            hint(2) = 0.5;
            
        % splitting by 2 protons
        elseif(npro==2)
            hint(1) = 0.25;
            hint(2) = 0.5;
            hint(3) = 0.25;
        
        % splitting by 3 protons   
        elseif(npro==3)
            hint(1) = 0.125;
            hint(2) = 0.375;
            hint(3) = 0.375;
            hint(4) = 0.125;
        
        % splitting by 4 protons
        elseif(npro==4)
            hint(1) = 0.0625;
            hint(2) = 0.25;
            hint(3) = 0.375;
            hint(4) = 0.25;
            hint(5) = 0.0625;
        
        % splitting by 5 protons
        elseif(npro==5)
            hint(1) = 0.03125;
            hint(2) = 0.15625;
            hint(3) = 0.3125;
            hint(4) = 0.3125;
            hint(5) = 0.15625;
            hint(6) = 0.03125;
        
        % Splitting by 6 protons
        elseif(npro==6)
            hint(1) = 0.015625;
            hint(2) = 0.09375;
            hint(3) = 0.234375;
            hint(4) = 0.3125;
            hint(5) = 0.234375;
            hint(6) = 0.09375;
            hint(7) = 0.015625;
        end
      end
      
return
end

